Plasmon excitation by charged particles interacting with metal surfaces 
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Recent experiments (R. A. Baragiola and C. A. Dukes, Phys. Rev. Lett. 76, 2547 (1996)) with 
slow ions incident at grazing angle on metal surfaces have shown that bulk plasmons are excited 
under conditions where the ions do not penetrate the surface, contrary to the usual statement that 
probes exterior to an electron gas do not couple to the bulk plasmon. We here use the quantized 
hydrodynamic model of the bounded electron gas to derive an explicit expression for the probability 
of bulk plasmon excitation by external charged particles moving parallel to the surface. Our results 
indicate that for each q (the surface plasmon wave vector) there exists a continuum of bulk plasmon 
excitations, which we also observe within the semi-classical infinite-barrier (SCIB) model of the 
surface. 



It is well known that charged particles interacting with 
solid surfaces can create electronic collective excitations 
in the solid. These are bulk 0] and surface [D plasmons. 
In the absence of electron-gas dispersion, the scalar elec- 
tric potential due to bulk plasmons vanishes outside the 
surface hence, in this case probes exterior to the solid 
can only generate surface excitations. That electron-gas 
dispersion allows external probes to interact with bulk 
plasmons was discussed by Barton [Q and Eguiluz [^, 
and more recently by Nazarov et al Q. Nevertheless, 
the fact that within a non-local description of screening 
bulk plasmons do give rise to a potential outside the solid 
has been ignored over the years , even when electron- 
gas dispersion has been included; also, current non-local 
theories of plasmon excitation by external charges have 
not shown evidence for bulk plasmon excitation outside 
the solid j^-p^ . Recently, Baragiola and Dukes have 
studied the emission spectra produced by slow ions that 
were incident at grazing angle; their data indicate that 
the bulk plasmon is importantly involved in the emis- 
sion process, though the projectiles are not expected to 
have penetrated into the solid. Bulk plasmon excitation 
in electron emission spectra produced by slow multiply 
charged ions has also been investigated with projec- 
tiles that may enter the solid. 

In this letter we derive, within the quantized hydrody- 
namic model of the bounded electron gas [p^ , an explicit 
expression for the probability of bulk plasmon excitation 
by external charged particles moving parallel to a jcUium 
surface. Our model, which assumes a sharp electron den- 
sity profile at the surface, neatly displays the role of bulk 
plasmon excitations in the interaction of charged parti- 
cles moving near a metal surface. We also demonstrate 
that our results for the total energy-loss probability agree 
with standard calculations [Q derived either by solv- 
ing the linearized Bloch hydrodynamic equations JTsI l or 
within the semi-classical infinite-barrier (SCIB) model of 



the surface p^Jl7[ | with the hydrodynamic approxima- 
tion for the bulk dielectric function. Though it has been 
generally believed that when charged particles move out- 
side the solid the energy loss predicted in these models 
is fully described by the excitation of surface plasmons 
§ , we demonstrate that for each q (the surface plasmon 
wave vector) excitation of both a discrete surface plas- 
mon and a continuum of bulk plasmons contribute to 
the total energy-loss probability, in agreement with the 
prediction of our quantized hydrodynamic scheme. 

Take an inhomogeneous electron system capable of self- 
oscillations about a ground state described by density- 
functional theory (DFT) In the hydrodynamic limit 
Jigf the system is characterized by the electron density 
n{r,t) and a velocity field u(r,t). The total energy of 
the system can then be expressed as (we use atomic 
units throughout, i.e., — Ti = rUe — 1) 



H [n(r, t)] + j dvn[r, t) 
I n{r, t)n{r' , t) 



drdr 



r - r' 



(1) 



where irrotational flow has been assumed, i.e., u(r,t) — 
—Vip, and retardation effects have been neglected. 
G[n{r,t)] represents the exchange, correlation and in- 
ternal kinetic energies of the electron system. Vq is the 
electrostatic potential due to the neutralizing ionic back- 
ground of density no, and Vi represents the external per- 
turbation. From Eq. (1) the basic hydrodynamic equa- 
tions can be derived, following Bloch's approach ||l9|, and 
they can be linearized in the deviation n — no from the 
equilibrium value to find the existence of self-sustaining 
normal modes of oscillation. 

We consider a charged particle moving with velocity v 
outside of a metallic surface, along a trajectory that is 
parallel to the surface, with a classical charge distribution 
given at r = (r||,z) by pcxt(r, i) = Zi5{y\\ - vt)S{z - zq). 



1 



We represent the ionic background by a jellium model 
(the jeUium occupying the half-space z < 0), and assume 
a sharp electron density profile at the surface. We also 
neglect exchange-correlation contributions to G[n\, and 
approximate it by the Thomas-Fermi functional [ p^ . In 
this approximation, for each value of q (the wave vector 
parallel to the surface) there exist both bulk and surface 
normal modes of oscillation with frequencies given by the 
following dispersion relations : 



B \2 



and 



S\2 



K) 



/3\2)l/2 



(2) 



(3) 



respectively, where ujp ~ (47mo)^/^ is the so-called 
plasma frequency and /3 represents the speed of propa- 
gation of hydrodynamic disturbances in the electron sys- 

/3/5gF, If being the Fermi 
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tem. We choose |21 
momentum. 

Now we follow Ref. to quantize, after linearization, 
the hamiltonian of Eq. (1) on the basis of the normal 
modes corresponding to Eqs. (2) and (3), which we shall 
refer after quantization as bulk and surface plasmons, 
respectively. We find 



H = He 



iff 



(4) 



where Hq represents the Thomas-Fermi ground state of 
the static unperturbed electron system [jl^ 
are free bulk and surface plasmon hamiltonians: 



Hff and 



^0 =^J2 [1/2 + alJt) a^^pit) 



and 



Ho' 



q,p>0 



-y 



1/2 + ^^^] &t(i) 



(5) 



(6) 



where and A represent the normalization volume and 
the normalization area of the surface, respectively, and 
where a^^^pit) and 6q(t) are Bose-Einstein operators that 
annihilate bulk and surface plasmons with wave vectors 
(q, p) and q, respectively. H^^' are contributions to 
the hamiltonian coming from the coupling between the 
external particle and bulk/surface plasmon fields: 



B/S 



drpcxt(r, i)0^/^(r, t) 



(7) 



(j)^/^{v,t) representing operators corresponding to the 
scalar electric potential due to bulk/surface plasmons. 
Outside the metal [z > 0), 



q,p>0 



and 



^(r,i) = -^E f'^') e'''"" ' (9) 



where f^^p{z) and fq{z) are bulk and surface coupling 
functions, respectively: 



[p4+p2(^2 +^2/^2) +^4/(4^4)] 



1/2 



(10) 



and 



W + 27?)] 



1/2 



(11) 



and where 7^ represents the so-called inverse decay length 
of surface plasmon charge fluctuations [||: 



1 r 



-I3q- 



(12) 



We derive the potential induced by the presence of the 
external perturbing charge as the expectation value of 
the total scalar potential operator E3]: 



y"'^(r,t) 



< ^o|(/.g(r,t) + 0g(r,^)|^o > 
<^'o|*o> 



(13) 



where 1^*0 > is the Heisenberg ground state of the in- 
teracting system and where 0||(r, t) and (l)ff{r, t) are the 
operators of Eqs. (8) and (9) in the Heisenberg picture. 
Our results j2^ reproduce previous calculations for the 
image potential [defined as half the induced potential at 
the position of the charged particle that creates it] of a 
static (v — 0) external charged particle which in the 
case of a non-dispersive electron gas (/? = 0) coincides 
with the classical image potential |2J] Vit^{z) — — (4z)^^. 

The energy loss per unit path length of a moving 
charged particle can be obtained as the retarding force 
that the polarization charge distribution in the electron 
gas exerts on the projectile itself ]25|] : 

- ? = - / drpcxt(r, t)Vt/™'^(r, t) ■ v. (14) 
ax V J 

By introducing here the induced potential, which we eval- 
uate from Eq. (13) up to first order in Zi, we find [p6| 



dE 
dx 



dq r da;^[P,^^+P,^J , (15) 



B / S 

where Pq,/j represent probabilities per unit time, unit 
wave number and unit frequency for the excitation of 
bulk/surface plasmons with wave number q and fre- 
quency Lo: 
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pB ^ r^2 _J_ 



dp , \ =-S[i^ - ^q,p) (16) 



and 



Try 9^ — 



(17) 



For comparison, we note that by solving the hnearized 
Bloch hydrodynamic equations the total probability per 
unit time, unit wave number and unit frequency for the 
external particle to transfer momentum q and energy w 
to the electron gas is found as: 



Im{~l/[c.2-2/3%(A, + g)]} _ 



where 



2gzo 



A„ 



(18) 



(19) 



and where (5 is a positive infinitesimal. The same result 
is obtained within either the specular reflexion (SR) 
| p6[ or the SCIB model of the surface, as long as 
the hydrodynamic approximation for the bulk dielectric 
function is considered. 

Within the various semiclassical approaches leading to 
Eq. (18) the role played by bulk and surface plasmons 
goes unnoticed. Furthermore, Eq. (18) shows no evi- 
dence of losses at the bulk plasmon frequency, and it has 
been generally believed that the energy loss predicted by 
this equation originates entirely in the excitation of sur- 
face plasmons. That this is not the case is clearly shown 
in Fig. 1. This figure shows -Pg,tj, as computed from Eq. 
(18) with = 15.8eV corresponding to the bulk plasma 
frequency of aluminum metal and with a finite damp- 
ing parameter 5 accounting for the finite lifetime of plas- 
mon fields. The parallel momentum transfer, the velocity 
and the distance of the particle trajectory above the sur- 
face have been taken to be q = 0.4 a.u., v = 2a.u. and 
zq = 1.0 a.u., respectively, and different values of 5 have 
been considered. One sees that loss occurs at the surface 
plasmon energy given by Eq. (3), while a continuum 
of bulk plasmon excitations occurs at energies (see 
Eq. (2)), which all are over uj^f^ = (w^ _^ ^2^2^1/2^ 
predicted by Eqs. (16) and (17). In the limit as 5 ^ 0+, 
both bulk and surface contributions to the total energy- 
loss probability of Eq. (18) exactly coincide with the 
predictions of Eqs. (16) and (17), thus demonstrating 
the full equivalence between our quantized hydrodynamic 
scheme and the more standard semiclassical approaches. 
Also, either introduction of both Eqs. (16) and (17) into 
Eq. (15) or replacement of the integrand in Eq. (15) by 
the energy-loss probability of Eq. (18) with (5 — > result 



in exactly the same total energy loss. As the damping 
parameter increases, surface plasmon excitation broad- 
ens to energies over w^g, and for 5 « Wp/10 bulk plasmon 
contributions to the energy loss go unnoticed. 

Figure 2 exhibits, by a solid line, the energy loss per 
unit path length versus the distance zq from the surface, 
as obtained from Eq. (15). In this case a proton [Zi = 1) 
moves with velocity v = 2 a.u. parallel to the surface of 
a bounded electron gas of density equivalent to that of 
aluminum {rg — 2.07) [ pTj . Separate contributions from 
the excitation of bulk and surface plasmons are repre- 
sented by dashed and dotted lines, as obtained with the 
use of Eqs. (16) and (17), respectively. One sees that 
for zq < 2 a.u. the contribution to the energy loss from 
the bulk channel is important, while for larger values of 
zq the surface channel alone gives a sufficiently accurate 
description of the total energy loss. 

Figure 3 shows, by a solid line, results for the energy 
loss of Eq. (15) , as a function of the velocity of the projec- 
tile, with Zi = 1, Ts = 2.07 and zq = 1 a.u.. Dashed and 
dotted lines represent separate contributions to the total 
energy loss from the excitation of bulk and surface plas- 
mons, respectively. We note that for a projectile moving 
at zq = 1 a.u. the contribution to the energy loss from the 
bulk channel is observable for all velocities, the relative 
importance of bulk plasmon excitations becoming more 
important at velocities around the plasmon threshold ve- 
locity when the projectile has enough energy to excite 
plasmons. For comparison, the result one obtains either 
from Eq. (17) or Eq. (18) in the case of a non-dispersive 
electron gas {(3 = 0), which coincides with the classical 
formula of Echenique and Pendry [ p8| , is represented by 
a dashed-dotted line. 

In conclusion, we have used the quantized hydrody- 
namic model of the bounded electron gas to demonstrate 
that bulk plasmons undergo real excitations, even in the 
case of charged particles that do not penetrate into the 
solid. We have derived explicit expressions for the prob- 
ability of both bulk and surface plasmon excitation by 
external charged particles moving parallel to the surface, 
which neatly display the role that bulk plasmon excita- 
tion plays in the interaction of charged particles mov- 
ing near a metal surface. The full equivalence between 
our quantized hydrodynamic scheme and the more stan- 
dard semiclassical approaches has been demonstrated. It 
has been generally believed that energy loss predicted by 
these approaches originates entirely in the excitation of 
surface plasmons. However, we have shown that for each 
value of the wave vector parallel to the surface both a dis- 
crete surface plasmon excitation and a continuum of bulk 
plasmon excitations contribute to the total energy-loss 
probability. We have also presented explicit calculations 
of the energy loss per unit path length of protons moving 
outside of a metal along a trajectory that is parallel to 
the surface, and our results indicate that the contribution 
from the bulk channel is important for all projectile ve- 
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locities, as long as the distance from the surface is smaher 
than a few atomic units. 
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FIG. 1. Energy-loss probability of Eq. (18), as a func- 
tion of uj/iiJp, with uip = 15.8eV, q = 0.4a.u., v = 2a. u., 
zo = la.u., Zi — 1 and various values of the damping pa- 
rameter 5: 10~*a;p (solid line), 10~^ci;p (dashed line), 10~^ujp 
(dotted line), and 10~^a'p (dashed-dotted line). 

FIG. 2. Energy loss per unit path length versus the dis- 
tance from the surface zo, as obtained from Eq. (15) with 
either the use of Eqs. (16) and (17) or the use of Eq. (18) and 
S —> (solid line). Separate contributions from Eqs. (16) and 
(17) are represented by dashed and dotted lines, respectively. 
A proton {Zi = 1) is assumed to move with velocity v = 2 a.u. 
parallel to the surface of a bounded electron gas with an elec- 
tron density equal to that of aluminum {rs — 2.07). The 
dashed-dotted line represents the total energy-loss obtained 
from Eq. (15) with no cut-off in the momentum integration 

{Qc OO). 

FIG. 3. As in Fig. 2, but now the energy-loss per unit path 
length of Eq. (15) is represented as a function of the velocity 
of the projectile for a given value of the distance from the 
surface: zq — la.u.. The dashed-dotted line here represents 
the result obtained in the case of a non-dispersive electron gas 
(/3 = 0). 
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